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Abstract 

This is a review of theoretical attempts to describe production of heavy quark 
^ , bound states in nucleus-nucleus collisions, in particular, the relative suppression of 

J> I J/ip and -0' production observed in these reactions. The review begins with a survey 

■ of experimental data for proton-induced reactions and their theoretical interpretation. 

Q . The evidence for additional suppression in nucleus-nucleus collisions is discussed and 

I various theoretical models of charmonium absorption by comoving matter are presented 

and analyzed. The review concludes with suggestions for future research that would 
. help clarify the implications of J/ip suppression in Pb + Pb collisions observed by the 

NA50 experiment. 



1 Introduction 



Ever since their discovery two decades ago, bound states of heavy quarks have served as one 
^ I of the primary probes of quantum chromodynamics. The spectroscopy of charmonium and 

bottomonium states provided crucial tests for the theory of forces among quarks, culminating 
in recent years in the precise determination of the strong coupling constant by means of 
lattice simulations. 

Because of their relatively small size, bound states of heavy quarks interact only rather 
weakly with other hadrons. As Matsui and Satz pointed out more than a decade ago, this 
property makes them excellent candidates as probes of superdense hadronic matter, created 
in collisions between nuclei at high energy |Q. As long as this matter were composed of 
hadrons, so they argued, it would only weakly affect simultaneously created heavy quark 
bound states {J/ip, ip' ■, T, etc.). On the other hand, if a color deconfined state, i.e. a quark- 
gluon plasma, were formed in the nuclear reaction the formation of heavy quark bound 
states would be severely suppressed, because the attractive color force between the quarks 
is screened by the plasma. Quantitative calculations of this effect at finite temperature 
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can be performed in the framework of lattice gauge theory. These calculations yield a 
state-specific critical temperature Tj^^ above which the bound state (z) dissociates into the 
plasma. Whereas the charmonium states dissociate soon above the critical temperature Tc 
of the QCD phase transition {T^^^^ ^ 1.2 T^), the bottomonium states survive up to higher 
temperatures {Tj^^^ ^ 2Tc). 

A systematic account of how these characteristic properties of heavy quark bound states 
may be put to use in the study of superdense matter, especially for the purpose of estab- 
lishing the existence of a new deconfined phase of QCD, was presented several years ago 
by Karsch and Satz [Q. Since then, new theoretical insight into the elementary process of 
heavy quarkonium production has led to a reassessment of the suppression of charmonium 
production in hadronic interactions with nuclei. In fact, as will be discussed in the next sec- 
tion, all experimental data taken with nuclear projectiles (^^O, ^^S) up to 1994 can probably 
be explained as final state interactions of the produced heavy quark pair with the nucleons 
contained in the colliding nuclei. The discovery of enhanced suppression of J/ip production 
in collisions between two Pb nuclei at the CERN-SPS in 1996, therefore, generated intense 
interest. Does the "anomalous" J/ip suppression observed in Pb + Pb collisions indicate the 
formation of a quark-gluon plasma in these reactions? Has this new state of matter been 
finally detected? 

In order to answer this question, all other conceivable mechanisms that could explain 
the observed suppression effect need to be ruled out. This is not an easy task because 
reactions among heavy nuclei at high energy are processes of great complexity, and there is 
little hope that they can be understood in all details. It is therefore necessary to allow for 
a considerable range of parameters in different reaction models and to exclude competing 
suppression mechanisms under a variety of assumed scenarios. 

It is the purpose of this review to assess the status of the theoretical tools that are 
needed to establish the possible origin of the observed "anomalous" J/ip suppression. The 
next section summarizes the current experimental information on J/tp and ip' suppression 
in proton-nucleus and nucleus-nucleus collisions. The explanation of all data (except those 
from Pb + Pb collisions) in terms of nuclear final state interactions will also be discussed 
in this section. In Section 3, a number of possible mechanisms for the enhanced suppression 
seen in Pb + Pb collision will be surveyed and their current theoretical uncertainties will 
be analyzed. Section 4 deals with the quantitative analysis of the Pb + Pb data in the 
framework of specific reaction models. The review closes with a list of suggestions for future 
research, mostly theoretical, that could help sharpen the tool of charmonium spectroscopy 
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for the exploration of the properties of superdense hadronic matter. 

The discussion here is mostly based on publications that were available in June 1997. 
Later developments are briefly discussed in a separate section at the end of the review. For 
a survey of the state of this research field from a somewhat different perspective, the reader 
is urged to also consult the recent review by Kharzeev 0. 

2 Phenomenology 

2.1 Proton-Induced Reactions 

The production of charmonium states in proton interactions with nuclei has been studied 

extensively at Fermilab as well as at CERN. Since these experiments observe the decay 

mode of the J/ip and ip', they also provide information on the nuclear target dependence 

of the fJ^^fi^ continuum spectrum produced by light quark- ant iquark annihilation (Drell- 

Yan process, DY). Whereas the Drell-Yan cross section is observed to grow almost exactly in 

proportion to the nuclear mass number A, the J/ip and ip' production cross sections grow less 

rapidly. The target dependence is often parametrized as a power law, a = ctqA"'. A value a < 

1 indicates "suppression" of charmonium production, compared with the naive expectation 

a/A = const, obtained when one neglects final state interactions. The argument is that the 

strictly linear A-dependence of the DY cross section rules out initial state interactions as the 

origin of the reduced production of charmonium states.^ 

The A-dependence of the nuclear J/ip and ip' cross sections is remarkably similar. Both 

can be fit by the same exponent a ~ 0.92, as shown in Fig. |I|. Since the mean radius of the 

ip' is much larger than that of the J/ip, the equally strong suppression indicates that the final 

state interactions occur with a state that is not an eigenstate of the (cc) system but rather 

with a common precursor of the J/tp and ip'. This does not come as a surprise, because the 

center of mass of the produced (cc) pair moves rapidly with respect to the target nucleus. 

As viewed from the rest frame of the (cc) pair, the nucleus is highly Lorentz contracted and 

thus the (cc) pair leaves the target nucleus before the components corresponding to different 

quantum mechanical eigenstates of the charmonium system have had time to decohere. 

^This argument is not strictly correct, because heavy quarks are predominantly formed by gluon fusion or 
gluon scattering. Gluons might be more susceptible to initial state interactions than quarks and antiquarks. 
There is, indeed, evidence for enhanced initial state scattering of gluons from the p^-spectra of heavy quark 
states produced in p+p reactions, giving rise to a broader "intrinsic" transverse momentum spread of gluons 

i- 
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Figure 1: Ratio of nuclear cross section 
compared to deuterium for Drell-Yan pairs 
(DY), J/ip and ip' production. Note that 
the nuclear suppression for J/ip and ip' is 
the same. Data from experiment E772 at 
Fermilab H. 



The influence of the nuclear target on the decoherence process has been studied ex- 
tensively in different theoretical approaches. However, the size of the nuclear suppression 
remained a great puzzle until quite recently. The problem is that the coherent superposition 
of eigenstates which is initially produced is characterized by a small geometric size of order 
~ 0.1 fm. On the other hand, the sublinear rise with A of the charmonium produc- 
tion cross sections requires a final state interaction cross section accN of the order of several 
millibarn. A comprehensive analysis of all available data on charmonium production in 
p + A collisions within the Glauber model yields a value for the absorption cross section of 
a^^J) = (7.3 ± 0.6) mb. 

2.2 The Color Octet Model 

The resolution of this apparent paradox can be achieved if one assumes that the (cc) pair is 
originally produced in a color octet state [0] which strongly interacts with its environment. 
The color octet model was originally motivated by unexpectedly large cross sections for 
charmonium production observed at the Tevatron. Perturbative calculations fell short of 
the measured cross section for J/tp and tp' production at high px by at least one order of 
magnitude. The observation ||^ that the data obtained by the CDF collaboration can be 
explained if high-p-p charmonia are predominantly formed by gluon fragmentation, led to the 
hypothesis that the charmonium wavefunctions contain a significant component where the 
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(cc) pair is in a color octet state accompanied by a soft (valence) gluon: 

mn'S^)) = 0(l)|ccpS«]) + Oiv,)\cc['PP]g) + Oivl)\cc['S[%g) + .... (1) 

Here Vc denotes the velocity of the bound charm quarks. This decomposition of the char- 
monium wavefunction can be unambiguously defined in the framework of the nonrelativistic 
expansion of QCD (NRQCD) developed by Lepage and collaborators 0. The higher Fock 
space components are suppressed by powers of the velocity vq of the heavy quark. In the 
heavy quark limit vq ~ ^asimq) -C 1, producing a well ordered expansion into components 
with an increasing number of valence gluons. The reason for the suppression of higher Fock 
space components is that "physical" gluons couple to the color current of the heavy quarks 
which is proportional to their velocity vq = pq/uiq. 

In the limit ^ 1, where the structure of the heavy quark bound states is essen- 
tially Coulombic, the momentum scale for the production of the heavy quark pair, 2mQ, 
is much larger than the momentum scale associated with the heavy quark bound state, 
Pqin^Si) = asTnQ/2n. In the case of the charmonium ground state, Pc ~ 1 GeV com- 
pared with 2mc ~ 3 GeV, which still provides for a reasonable separation of scales. The 
total S-matrix element for the production of a cc[n^S'i] state may therefore be factorized 
into a perturbative matrix element for the elementary production process, such as gluon 
fusion {gg —>■ cc) or gluon fragmentation {g* cc), and a nonperturbative matrix element 
Ai{cc — i> ipi'i^^Si)) describing the evolution of the elementary (cc) pair into the strong inter- 
action eigenstate (1), as illustrated in Fig. |^. There are different matrix elements, Aii and 
A^8, depending on the color channel in which the (cc) pair is produced. 

Figure 2: Schematic representation of the 
color-octet formation mechanism for the 
charmonium states J/ip and ip'. The inter- 
mediate state ^S[^^ is a color octet which 
combines with a soft gluon to form the 
color-singlet charmonium state. 

Assuming that charmonium production at high p^p at the Tevatron is dominated by the 
color octet channel, the matrix elements Ais for the formation of various eigenstates of the 
charmonium system can be deduced empirically [|1^, 0]. Their values are found to scale 
in accordance with the counting rules, the powers of Vc, predicted by (1). Although the 
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color octet model yields an excellent description of the measured differential cross sections, 
one has to caution that crucial predictions, such as the strong transverse polarization of 
J/%1) produced at high p^, have not yet been confirmed. It should also be noted that the 
confirmation of this prediction would not automatically prove that the color octet mechanism 
also dominates J/%1) production at small pt, where the cross section is concentrated in nuclear 
collisions. 

An important theoretical problem is how quickly the color octet state neutralizes its 
color charge. Since the neutralization mechanism is, by assumption, a nonperturbative soft 
process, it is not rigorously calculable, but some qualitative estimates are possible. Since the 
force between the (cc) pair in a color octet state is repulsive, the formation of a bound state 
requires that the relative (cc) wavefunction must be a color singlet at distances comparable 
to the average radius of the charmonium state. More precisely, the relative magnitude of the 
octet and singlet components as a function of separation r should behave as 



A crude estimate is obtained by replacing \Ecc — Vg] with the binding energy of the charmo- 
nium state, Eb{v?Si), yielding the neutralization time 



Clearly, a more quantitative estimate within the framework of a coupled-channel potential 
model for charmonium would be very useful. In the laboratory frame, the octet-to-singlet 
transition time is Lorentz dilated due to the rapid motion of the (cc) pair with respect to the 
target nucleus. Since the average traversed length of nuclear matter is of the order of 5 fm for 
a heavy nuclear target, the precise value of rg^i becomes important for charmonium states 
produced at midrapidity. However, in first approximation, it is reasonable to expect that the 
color neutralization occurs predominantly outside the target nucleus in p -\- A experiments 
performed at several hundred GeV/c beam energy. 

2.3 Nucleus-Nucleus Collisions 

The production of J/ip and ip' in nucleus- nucleus collisions has been studied in the experi- 
ments NA38 and NA50 at CERN. NA38 investigated p + W, p + U, ^^O + U, and ^'^S + U 
collisions at 200 GeV/nucleon; NA50 took data for ^ospb -\- Pb coUisions at 158 GeV/nucleon. 






(3) 
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Figure 3: Effective nuclear absorption cross 
section in the Glauber approximation for 
J/ip deduced from the NA38 data. The 
heavy ion data (black squares) fall on the 
same trend as the data from proton induced 
reactions. 



Figure 4: Illustration of the definition of L, 
the equivalent thickness of nuclear matter 
sweeping over the location of the created 
(cc) pair. 



The results of experiment NA38 showed that the production of J/ip in ^^O- and ^^S-induced 
reactions follows the nuclear mass dependence expected fromp + A collisions see Figure (|^). 
The total production cross section ajj^ grows like0 

aj/^ = a^]jl ■ A,A, exp (-aJSV^^l^i, ^2)) = a^^JflA^A^S^, (4) 

where L{Ai,A2) is the average length of nuclear matter traversed by the (cc) pair after its 
formation at some moment during the collision. The definition of L is illustrated in Fig. ^ 
In the framework of Glauber theory, the nuclear suppression factor Sn (and hence 
L{Ai,A2)) is defined as 

Sn = j d% (fh'dzdz'pA, {h\ z)pa, (b - b', z')Ta, {z, b')TA, {z\ b - b') (5) 

with the nuclear profile function 

Ta{z, b) = exp (-(A - l)aip dz'pA{b, z')^ . (6) 

^Note that p^L is a Lorentz invariant expression, which counts the area density of nucleons interacting 
with the (cc) pair after its formation. 
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Here pA{b, z) is the nuclear density a'^f^'' distribution, normalized to unity so that Sn = 1 
for p + p collisions. 

The fact that nuclear production of J/ ip appears to follow the Glauber model prediction 
(^) may be interpreted to mean that nothing "abnormal" occurs in nuclear reactions up to 
projectile mass Ai = 32. This argument was first presented in this form by Capella, et al. 
and by Gerschel and Hiifner [^. It is worth noting that their conclusion originally 



was quite controversial because the absorption cross section a[^^^ was expected to be much 
smaller than the required 6-7 mb. This state of affairs has been radically changed by the 
introduction of the color-octet model that provides a natural explanation for cr^g^'*. 

The results obtained for ip' production on the other hand, cannot be explained by (Q). 
The ip' cross section is found to be more suppressed in ^^S + U reactions than the J/tp cross 
section. This is in marked contrast to p+A collisions where the ratio cr^/ / ctj/^ is independent 
of the target mass. The different behavior of the ip' becomes even more evident when the 
production of charmonium states is studied as a function of the total transverse energy Ej' 
carried by the collision fragments. Within uncertainties, is a measure of the impact 
parameter of the nuclear reaction. Selected ranges of E^, therefore, correspond to different 
values of the average target thickness L{Ai, A2) in (^. The relation between L and E^ can 
be modeled in the Glauber approximation or some other geometric collision model. The 
relative suppression of ip' production compared with J/ ip production increases with growing 
L (or Et)-, as shown in Fig. ^ . On the other hand, (Jj/^ in different i^^-windows follows 
exactly the Glauber model formula (^), if the i?T-dependence of L is taken into account. 

The stronger suppression of the ip' indicates the presence of a new suppression mechanism 
in S + U collisions that is absent m p + A collisions. One candidate for this mechanism is 
absorption by "comovers", i.e. by secondary hadrons produced at about the same rapidity 
as the ip' . Since they are slowly moving with respect to the ip' , there is no time dilation of 
the formation time of the ip' for interactions with these secondary particles. We will return 
to this issue in section 4, when we discuss models for comover absorption. 

In Pb + Pb collisions at 158 GeV/nucleon, as shown in Fig. ^, also the cross section 
for J/ip production is found to be more strongly suppressed than predicted by The 
amount of additional suppression increases with L (or Et)-, similar to the effect observed 
for ip' production in S + U collisions. For ip' production in the Pb + Pb system, strong 
additional suppression is observed even in the lowest Et range, as shown in Fig. |^. The 
suppression factor for the highest Et (or L) window is only insignificantly larger than that 
found in the most central S + U collisions. Taken together, the S + U and Pb + Pb data 
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Figure 5: The cross section for ^0' produc- 
tion is suppressed relative to that for J/ip 
production in S + U coUisions. The sup- 
pression increases with L. Heavy ion data 
from experiment NA38. 



indicate the onset of a new, "abnormal" suppression mechanism for charmonium production, 
first for the ip' and later for the J/ip, which appears to approach saturation for the ■0' in the 
most central Pb + Pb events. 

The most obvious difference between p + A, S + U, and Pb + Pb interactions at SPS 
beam energies is the number of secondary particles. As already mentioned before, a sig- 
nificant fraction of the secondary hadrons are produced in the same region of rapidity as 
the charmonium. Is it possible to explain the "abnormal" suppression as absorption of the 
J/ip and ■0' on these comoving hadrons? At first glance, this explanation seems eminently 
plausible. The absorption cross section of the ip' and comovers should be much larger than 
that of the J/t/j, both because of the larger geometrical size of the ip' and its lower energy 
threshold for dissociation. Therefore, the comover effect should set in much earlier for the ■0' 
than for the J/ip, just as observed. The question is thus decidedly of a quantitative, rather 
than qualitative, nature: 

1. Is it possible to explain the A- and -Ey-dependence of the "anomolous" part of J/ijj 
and '0'-suppression by a common set of parameters? 

2. Do the required parameters, i.e. comover absorption cross sections, agree with theo- 
retical estimates or values deduced from other relevant data? 

The following sections address these two questions. We shall first review what theory 
can tell us about cross sections for the absorption of charmonium on hadrons made of light 
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Figure 6: Results from NA50 for the sup- 
pression of J/ production in Pb + Pb col- 
lisions as function of the Glauber parameter 
L. (From 0.) 



Figure 7: Results from NA50 for the sup- 
pression of production in Pb + Pb colli- 
sions as function of the Glauber parameter 
L. (From 0.) 



quarks. We shall then turn to phenomenological models that aim at finding a consistent set 
of parameters describing the S -|- U and Pb + Pb data. 

3 Charmonium Absorption by Comovers 
3.1 Overview of Mechanisms 

Most mechanisms that have been proposed to describe the absorption of J /if) or ip' on 
comoving hadronic matter fall into three categories:^ 

1. Deconfinement: If the comoving matter density is sufficiently high, the matter may be 
in the form of a quark-gluon plasma. In this case the charmonium states cannot form 
at all, as first argued by Matsui and Satz, because the attractive color force between 
the (cc) pair is screened by the plasma. The color-singlet (cc) potential then has the 
form 

Kc(r) = (7) 
6 r 

■^A comprehensive description of QCD based theoretical approaches to J/ip interactions in hadronic matter 
can be found in ||l^ . 
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where both as and are functions of the density of the plasma. For a thermahzed 
plasma with no net baryon excess, finite temperature perturbation theory predicts 

/i^ = (l + ^iVf) g'T' (8) 

where A^f = 3 is the number of light quark flavors. Monte-Carlo simulations of the 
SU(3) lattice gauge theory at finite temperature |TB|, show that the inverse color 
screening length /i in the range < T < 3Tc can be parametrized as ~ 2.2T. For 
every bound state of a heavy quark- ant iquark pair then exists a critical temperature 
Td, above which the state disappears. Detailed studies p| have revealed that (■?/'') ^ 
Td(Xc) ~ Te, Td(J/V^) ^ l.2T„ and Td(T) = 2T,. 

Thermal dissociation: If the (cc) pair is immersed in a thermal hadronic gas, its internal 
modes will eventually also become thermally occupied. All those states above the DD 
threshold will dissociate. As the temperature rises, the thermal fraction of dissociated 
states increases until, eventually, the bound states make up a negligible contribution. 
In condensed matter physics, the temperature above which almost no bound states 
remain is called the Mott temperature. In this picture, the J/ ip and ip' could essentially 
disappear even if quarks and gluons remain confined at any temperature. However, 
the critical question is how fast the internal states of a heavy quark pair become 
thermalized by interactions with a hadronic medium. This is an issue of kinetics, 
rather than thermodynamics, and will be addressed in the next subsection. 

Pre-thermal dissociation: In this scenario, charmonium states could be excited above 
the dissociation (DD) threshold by interactions with a medium of comovers that had 
no time to equilibrate. Two popular examples are (i) a pre-thermal parton plasma 



Ig, (ii) a system of color flux tubes [2^, |21[| . 

The first scenario, a pre-equilibrium parton plasma, obviously only makes sense if the 
dissociation process is faster than the thermalization of the plasma: r^is < Tgq.. This 
seems an unlikely case, because parton thermalization is governed by the cross section 
(jgg. But in the perturbative realm of QCD, cigg ~ |o"gc, hence the plasma should 
equilibrate faster than a (cc) state can be dissociated. 

Nonetheless, quantitative model studies of the competition between equilibration of 
partons and J/ ip dissociation would be useful. These could be carried out at two levels. 
One could study the dissociation of a fully formed charmonium state by gluon impact in 
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the environment of a gluon distribution that itself is equihbrating. Or one could study 
the influence of a prethermal gluon bath on the conversion process from color octet 
(cc) to the color singlet charmonium state. The latter would require some modeling 
of the nonperturbative matrix element Aig, but a semi-quantitative treatment of this 
mechanism should be possible. 

The second scenario, dissociation by color flux tubes, again requires that the decay rate 
of flux tubes by light quark pair production, Tgg, is slower than that of the dissociation 
rate of the charmonium state. This also appears unlikely, because QCD flux tubes 
are known to break rapidly (r^g < 1 fm/c), while a (cc) pair cannot be pulled apart 
too rapidly by virtue of the large mass of the c-quark. Again, simple model studies 
of the competition between these two processes would be easily carried out, and the 
dependence on heavy quark mass investigated. In this scenario, too, one can study the 
possible influence of the presence of coherent gluonic fields on the color octet-to-singlet 
conversion process. 



3.2 Thermal Dissociation Kinetics 

Since, in general, the cross sections cJhh among light hadrons are much larger than those 
between light hadrons and charmonium, a^^, it also makes sense to restrict the study of 
dissociation of charmonium by light hadrons to a thermal hadronic environment .0 As we 
argued before, the same consideration applies to a partonic environment above Tc. 

There are two possible approaches to the dissociation problem: one at the quark level 
where the large mass of the c-quark is used to separate perturbative from nonperturbative 
aspects of the problem, and another one that makes use of effective hadronic interactions. In 
this section, we will consider the constituent (cc)-approach. An effective hadronic Lagrangian 
approach will be discussed in the following section. 

As first explored in detail by Peskin and Bhanot and Peskin interactions between 



heavy quark bound states and light hadrons can be described perturbatively, if the heavy 
quark mass itlq is large enough. The reason is that the characteristic momentum scale of 
the heavy quark bound state is of order ctg^q which is large. As a result, only gluons with 
momenta larger than agrriQ can resolve the internal color structure of the {QQ) bound state 
and interact with it. To gluons with much smaller momenta, the (QQ) state appears as an 
inert color singlet. The small size of the (QQ) state then allows for a systematic multipole 
*The validity of this assumption can be checked with the help of hadronic cascade models. 
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expansion of its interaction with external glue fields, where the dipole interaction dominates 
at long range. 

The potential problem with this approach is that the c-quark is not quite heavy enough. 
For instance, the binding energy of the J/ip is about 650 MeV, which is neither small com- 
pared to the reduced mass \mc ~ 700 MeV of the charmonium system, nor large on the 
scale of light hadron masses. Gluons with momenta of order 1 GeV/c or even slightly less 
can effectively resolve the internal color structure of the J/V', but it is not clear that they 
belong in the perturbative domain of QCD. 
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Figure 8: Forward scattering between a light hadron and & J /ip is, in leading order, described 
as two-gluon exchange. The total absorption cross section is related to the imaginary part 
of this forward scattering amplitude. 

By virtue of the optical theorem, the absorptive cross section between a light hadron h 
and a charmonium J/ip can be expressed as the imaginary part of the forward scattering 
amplitude where the interaction is dominated at long distances by two-gluon exchange (see 
Figure |^). In this picture, light hadrons interact with the J/ip only via their glue content. 
Kharzeev and Satz ||2^ applied the Bhanot-Peskin formalism to tt — J/ip scattering. Using 



the operator product expansion to separate long-distance physics from the (assumed) pertur- 
batively calculable structure of the (cc) system, they found that the absorptive cross section 
is proportional to the part of the gluon structure function Gh{x) of the interacting hadron 
that is sufficiently energetic to dissociate the J/ip. Since the gluon structure functions of 
light hadrons typically are soft, i.e. Gh{x) (1 — x)^, only highly energetic hadrons are 
capable of exciting a. J/ip above the dissociation threshold. The peak cross section^ for J/ ip- 
dissociation by a gluon in the dipole approximation is about 3 mb at a gluon momentum 
around 1 GeV/c. Due to the softness of the gluon structure function, the pion momentum 
must reach 5 GeV/c before attaining an absorption cross section in excess of 1 mb. The 

^The fact that <j''g^^'^ peaks and then falls off is presumably an artifact of the dipole approximation. If 
higher orders were retained, the cross section should approach a constant value at high energy, 2a gc- 
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thermally averaged cross section, in this framework, remains less than 0.1 mb for a pion gas 
within any realistic temperature range. 




Figure 9: Inelastic scattering of a light Figure 10: The J/V" state is dissociated into 

hadron on a J/il'-, yielding two D-mesons a (cc) octet pair by absorption of an ener- 

in the final state. getic gluon from its environment. 



The calculation of pj] is clearly incomplete, because pions are not very efficient at disso- 



ciating the J/ip. Heavier hadrons, such as the p-meson, can make significant contributions 



although they are much less abundant than pions. Moreover, the calculation |Q does not 
account for the energy balance between initial and final states, i.e. that the total mass of the 
colliding hadron is available to convert the J /if) state into a pair of D-mesons (see Fig. |^). 
This effect is especially important for heavier mesons, such as the p. 



Repeating the calculation ||2^ in a somewhat more general form, one finds that the 
dissociation rate of J/T/'-mesons in a (thermal) hadronic medium is given by (see Appendix 
A for details) 

1 

Tdis^-TT du{\Eico)\^)7ra,a^, (9) 

where a is the Bohr radius of the J/ip and -^{\E{uj)\'^) denotes the time-average fiuctuation 
density of color (electric) fields with frequency uj in the hadronic environment, as illustrated 
in Fig. |Ty. The lower limit of the integral accounts for the dissociation energy threshold, and 
the factor (ttsO^) arises from the coupling to the color dipole moment of the (cc) pair in the 
J/ip. Equation (|^) is applicable to any kind of hadronic medium, whether thermal or not. 
For example, one could estimate the color field fiuctuations in the framework of a random 



field model as proposed by Hiifner et al. [^. For a thermalized medium, QCD sum rules 



or lattice gauge theory could be applied. Because some heavy quarkonium states, especially 
the T do not disappear right at the critical temperature, the color field fiuctuation spectral 
density is also an important quantity to know in the deconfined phase, where it could be 
evaluated either in perturbation theory or by lattice simulations. 
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Figure 11: In the framework of the constituent 
quark model, the absorption of a J/?/^ on a hght 
quark meson is viewed as quark exchange reaction. 
At lowest order, four one-gluon exchange diagrams 
contribute to the amplitude. (From Martins et al. 



The real limitations of lie in the fact that the J/ ip is not truly a perturbative bound 
state but probes also the confining part of the cc-potential. It is therefore of interest to 
investigate the dissociation rate also in the framework of other approaches that model quark 



confinement. One such model is the constituent quark model of Isgur and Karl p6[. In this 
model the reaction 

n + J/tP ^ D* + D {D + D*) (10) 

is viewed as quark exchange where a c-quark and a light quark change sides. The cross 
section for this reaction was calculated by Martins, Blaschke, and Quack in the first Born 
approximation (see the diagrammatic representation in Fig. |ll]). Including D*D, DD*, and 
D*D* final states, the cross section peaks around 1 GeV kinetic energy (cm.) at a value 



of 15 mb (see Fig. |TT]). The magnitude of this cross section is solely due to the action of 
the confining interaction between the quarks, which is modeled as a "color-blind" attractive 
interaction between the quarks which has a Gaussian momentum dependence. Because this 
interaction is taken as attractive independent of the color quantum numbers of the affected 
quark pair, it does not cancel for the interaction of a light quark with a point-like (cc) 
pair in a color-singlet state, in contrast to the one-gluon exchange interaction. Moreover, 
the magnitude of the obtained cross section invalidates the use of the Born approximation, 
unless a number of different partial waves contribute. This is impossible to ascertain from 
p7| , because no partial wave decomposition of the dissociation cross section is presented 
there. 

Since the results obtained within the constituent quark model differ by orders of mag- 
nitude from those obtained in the Bhanot-Peskin approach, it is useful to consider the 
absorption process also in an entirely different framework. This will be done in the next 
subsection. 
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Figure 12: Energy dependence of 
the J/ absorption on pions for the 
final states DD* and D*D*. (From 
the preprint version of Martins, et 
al. 



Note that the figure in the 
printed version differs shghtly from 
the one reproduced here.) 



3.3 Hadronic Dissociation: Effective Theory 

The most abundant mesons in a hot hadronic gas are tt, K, and p. These can induce the 
following dissociation processes when encountering & J/ip particle: 

Ti + J/i) ^ D + D\ D + D* (11) 
K + J/ij^Ds + D*, Ds + D\ D + Dl, D + (12) 
p + J/tp ^ D + D, D* + D* (13) 

The kinematic thresholds for these reactions are 

(Q) 640 MeV, (0) 385 MeV, (|I3|) -135 MeV (+155 MeV). 

The reaction p + J/tp D + D has the lowest total invariant mass threshold and is, in fact, 
exothermic. At thermal equilibrium it is only suppressed due to the relatively high mass of 
the p-meson which causes p-mesons to be less abundant than pions in a thermal hadron gas. 

From a microscopic point of view, all reactions listed above can be understood as quark 
exchanges, where the J/ip transmits a charm quark to the light meson and picks up a 
light (m, d, or s) quark. Since similar reactions among light hadrons typically occur with 



large cross sections, one can expect that the reactions (|TT]JT^,[T3|) also proceed with significant 
strengths above their respective kinematic thresholds. This is, indeed, the result obtained by 
Martins, Blaschke, and Quack which was discussed at the end of the previous subsection. 



Here, we follow a different approach |2^. In the effective meson theory, the exchange of 
a (eg) or (eg) pair, where q, q stands for any light quark, can be described as the exchange of 
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a Ds, D*, or D* meson between the J/ip and the incident hght meson (see Fig. jTS]). Near 
the kinematic threshold, where only a small number of final states are open, the description 
in terms of an effective meson exchange is expected to be valid. 
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Figure 13: Lowest order Feynman diagrams contributing to the charm exchange reactions 
(p!T|) and (p!3|) . Single lines represent pseudoscalar mesons; double lines denote vector mesons. 



In order to calculate the various Feynman diagrams for the reactions (|TT|,0,0) we need 
to construct the effective three-meson vertices. We do so by invoking a strongly broken U(4) 
flavor symmetry with the vector mesons playing the role of quasi-gauge bosons. Denoting 
the 16-plet of pseudoscalar mesons {'7i,ri,ri' , K, D, Ds,ric) by $ = (f)iTi, where the Tj are the 
f/(4) generators, and the vector meson 16-plet {u, p, (j), K* , D* , D*,ip) by = Vf'Ti, the 
free meson Lagrangian reads 

£o = tr(9'^<l>t 5^$) - tT{d''V^'''){d^V^ - d^V^) - tr(<l>tMp<l>) + ^tr(V'^tMvV^). (14) 

Here Mp and My denote the mass matrices for the pseudoscalar and vector mesons, respec- 
tively. Because of the heavy mass of the charm quark, Mp and My break the U(4) symmetry 
strongly down to U(3), the mass of the strange quark introduces a weaker breaking to U(2), 
and the axial anomaly further breaks the U(2) symmetry to SU(2) in the case of the pseu- 
doscalar mesons. All these symmetry breakings are embodied in the physical mass matrices. 
It is convenient, in the following, to work with the mass eigenstates. 

The meson couplings are obtained by replacing the space-time derivatives by the 
"gauge covariant" derivatives 



(15) 
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In first order in the coupling constant g this procedure leads to the following interactions 

Ant = ig ti {<^^V^^d^<^ -df'<^^V^<t>) 

+ig tr (9^Vt^[V^, VA - [V^^ V^I^^V.) . (16) 

If the U(4) flavor symmetry were exact, we would expect all couplings given by the same 
constant g. In view of the significant breaking of the flavor symmetry we anticipate that the 
effective coupling constants for different 3-meson vertices will have different values. We will 
see below to what extent this is true. 

In order to describe vr- and p- induced J/ip dissociation, we need the following vertices: 
ipDD, ipD*D*, TiDD*, pDD, and pD*D*. From (|l^) we derive the following interactions: 



^pDD = ^9pDDPi (pTid^D - d^DTiD) 
CpD'D' = -'-gpD*D*pf{D*''r,d,D:-id,D*nnD:). (17) 

The coupling constants g^oD, 9^)0*0*, QpDD and gpD*D* can be derived from the D and D* 
electric form factors in the framework of the vector meson dominance (VMD) model. 

If 7v denotes the photon-vector meson V mixing amplitude, the standard VMD analysis 
yields the relations 

2 

Ipfp = eml, 7^/^ = -eml (18) 
where universality of the vector meson couplings is assumed: 



fp — dpDD — gpD'D*, fip — gi)DD — g^D*D*- (19) 

The photon mixing amplitudes 7v can be determined from the leptonic vector meson decay 
widths [g9[: 

1 7v 

Tvee = -a^. (20) 
6 my 

Inserting the experimental numbers, we find 

fp ~ 5.6, U ~ 7.7. (21) 
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As is well known, different ways of deriving the value of the "universal" p-meson coupling 



fp yield values differing by about 20%. The coupling constants (pT]) must therefore be 
considered to be given with an error of this order of magnitude. 

Finally, the vr-meson coupling between D and D* can be obtained from the decay width 
of the D*-meson: D* Dn. Unfortunately, only an upper bound for this decay rate is 
known at present, corresponding to q-kdd* < 15. Theoretical estimates for this decay rate, 
based on QCD sum rules yield 



g^DD* ~ 8.8. (22) 

We adopt this value here. 

At the hadronic level, the charm exchange reaction between the J / ip and a light hadron 
can proceed either by exchange of a or a D*-meson. Here we shall only consider the 
D-exchange reactions, because they have an acceptable high-energy behavior even in the 
absence of form factors for the vertices. The D*-exchange reactions have cross sections that 
rise rapidly with energy due to the exchange of longitudinally polarized D*-mesons. This 
growth characteristic of the exchange of massive vector bosons is obviously unphysical and 
will be cut off by the mesonic form factors long before it becomes significant. However, 
this would require that we introduce an unknown parameter (the characteristic momentum 
scale of the internal mesonic structure). In order to avoid this complication, we here neglect 
D*-exchange cross sections. One final remark on D*-exchange: Regge theory dictates that, 
in the high energy limit, the charm exchange reaction is dominated by the exchange of the 
D*-trajectory. However, here we are not interested in charm exchange at high energies but 
near the kinematical threshold, because the relative motion of the hadrons is limited to 
thermal momenta. 

Analytical expressions for the amplitudes and cross sections for the processes correspond- 
ing to the Feynman diagrams in Fig. |13| are given in Appendix B. The charm exchange cross 
sections are plotted in Fig. as functions of the center-of-mass energy y/s. The + J/ip 
cross section (dashed line) starts at zero, because the reaction is endothermic, whereas the 
p + J /ip cross section (solid line) is finite at the threshold because of its exothermic nature. 
Of course, the reaction rate vanishes at threshold in both cases. Note that the cross sec- 
tions for these two J / ip absorption reactions are of similar magnitude over the energy range 
relevant to a thermal meson environment. 



Figure IT5| shows the pion and p-meson densities in an ideal thermal meson gas as a 



function of temperature. In the temperature range where a hadronic gas is likely to exist. 
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Figure 14: Cross sections for the charm ex- 
change reactions described by the diagrams 



of Figure 0, as functions of cm. energy. 
Dashed hne: pions, sohd hne: p-mesons. 



Figure 15: Thermal pion and p- meson den- 
sities in an ideal hadron g clS clS function 
of temperature. Dashed line: pions, solid 
line: p-mesons. The dotted line represents 
the density of pions above the dissociation 
threshold for J/ip. 



T < 170 MeV, pions are far more abundant than p-mesons, but most of these pions do not 
have sufficient energy to initiate the dissociation of a J/ip- In fact, if one counts only those 
pions and p-mesons above the kinematical threshold for J/?/'-dissociation, vr- and p-mesons 
are about equally abundant, or rather rare, because the effective density does not exceed 0.1 
fm~^ even at T = 200 MeV, where the hadron picture of a thermal environment probably 
already fails. Figure |1^ shows the J/ ip absorption rates in a thermal meson gas, as a function 
of temperature. Even at the (unrealistically) high temperature T = 300 MeV, the thermal 
dissociation rate is still so small that it corresponds to a lifetime around 10 fm/c. Thermal 
dissociation at T < 200 MeV is completely negligible on the time scale of the lifetime of a 
hot hadronic gas state in nuclear collisions. 

One can ask the question whether a form factor should be included in the meson vertices 
used to evaluate the Feynman diagrams in Fig. |T^. In principle, the concept of vector meson 
dominance assumes that the off-shell behavior of the vector meson propagators describes the 
form factor of the pseudoscalar mesons correctly, when the pseudoscalar mesons are on-shell. 
In the diagrams of Fig. O, however, the exchanged D-meson is off-shell by a considerable 
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Figure 16: Thermal J/ip absorption rates 
as function of temperature. The pion and 
p-meson rates are shown separately. 



Figure 17: Thermally averaged J absorp- 
tion cross sections as function of tempera- 
ture. The pion and p-meson rates are shown 
separately. 



amount, whereas the vector mesons all appear on-shell either in the initial or final state. It 
is unclear whether an additional form factor is needed in this situation. In any case, our 
result represents an upper limit, because an additional form factor would reduce the charm 
exchange rates. (E.g. a Gaussian form factor exp(— Q^/Qq) with Qo = 1-5 GeV would 
further reduce the absorption rates by slightly more than one order of magnitude.) 

As we noted above, our analysis also remains incomplete because of the neglect of D*- 
exchange reactions. It would be most interesting to consider also these reactions in the 
framework of a complete effective Lagrangian describing the interactions of pseudoscalar 
and vector mesons with quark content (gg), (Qg), (gQ), and (QQ), where g stands for 
any light quark flavor and Q denotes a heavy quark. Such a Lagrangian embodying chiral 
symmetry for the light quarks has recently been given by Chan . Loop corrections to the 
tree diagrams would also permit the study of form factor effects in this approach. 



4 Models versus the Data 
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4.1 Glauber Theory 

Probably the most conservative approach to calculate nuclear effects of suppression of char- 
monium production is based on the general framework of Glauber theory. This was, in fact, 
the starting point of the early analysis of Gerschel and Hiifner [|13] , and it forms the basis of 



more recent analyses of the available p + A and A + A data & 32 1. In the Glauber formal- 



ism the suppression factor Sn due to absorption on nucleons is given by @. The effect of 
absorption by comoving secondary hadrons is included by multiplying the integrand in (§) 
by an additional exponential absorption factor: 

o = I a u C 

with 



I ''''''' ^^-^'^''^ ^''^ 
-f dN,,{b,b') 



Tco{b,b') = exp 



-a5 



dy 



-dr 



(24) 



Here dNco{b, b')/dy is the comover density per unit of rapidity at impact parameter b, b' with 
respect to the two nuclear centers, and cr^f = is the effective J/ifj absorption cross 

section by comovers. The effect of quark deconfinement (complete absorption) can be easily 
described in this framework by setting a^^ = oo for dN^o/dy > dN^^^^ /dy. Tq corresponds to 
the formation time of the dense comover gas (typically taken as Tq = 1 — 2 fm/c), and 
denotes the time of hadronic freeze-out. 

It is reasonable to assume that the comover density is proportional to the transverse 
energy dEx/dy produced in a nucleus- nucleus collision at impact parameter b. A detailed 
analysis of dEx/dy measurements in S -|- U and Pb + Pb collisions shows that they can 
be well described by assuming that the transverse energy production is proportional to the 
number of participant projectile and target nucleons, as calculated in a geometrical picture 
or within the Glauber model itself. This is in agreement with a large body of data indicating 
that the "wounded nucleon" model ||3^ provides a good description of global observables in 
nuclear collisions. Using Bjorken's formula for the comoving energy density p5 



where vri?^ denotes the effective nuclear geometric cross section, and setting tq = 1 fm/c, 
the initial energy densities for different impact parameters in Pb + Pb and S + U collisions 
are found as shown in Fig. [l^. Clearly, for most impact parameters the density achieved 
in Pb + Pb collisions exceeds that produced in even the most central S + U collisions. 
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Figure 18: Initial energy densities according to the Bjorken formula (^) are shown as 
functions of impact parameter for (left) S + U and (right) Pb + Pb colhsions. (From 



ref. |6|. 



This observation opens the door |36| for an explanation of the anomalous J /if) production 
observed in Pb + Pb colhsions as the effect of deconfinement if a critical energy density 
ec ~ 2.75 GeV/fm'^ is exceeded (more precisely, if eoTo exceeds 2.75 GeV/fm^c). 

The earliest systematical analysis of suppression by comovers is due to Gavin and Vogt 
||37| . Using the parameters cr^g^^ = 4.8 mb, 0"°^ = 3.2 mb, and tq = 2 fm/c that were 
originally derived from the S + U data they find a good overall agreement with the 
measured suppression in Pb + Pb collisions P3. However, it must be noted that, according 



to our present understanding, their value of cr^g^^ may be too small to consistently explain 
the p + A data [Q, and is most likely much too large. The point here is that at tq = 2 
fm/c the formation of a color singlet (cc) state should be completed and, hence, the co mover 
absorption should be of hadronic size (ct^q < 1 mb) as discussed in the previous section. 

On the basis of their Glauber model analysis oip + A data on J/ ip suppression, Kharzeev 
et al. argue that there is no room for absorption by comovers in S + U collisions (left 
part of Fig. |19D. In order to explain the significant additional suppression observed in the 
Pb + Pb data (right part of Fig. |l^), these authors must invoke a suppression mechanism 
that sets in abruptly and strongly if dEx/dy, or rather the area density of wounded nucleons, 
exceeds a certain critical value. 



This is also found in the analysis of the same data by Armesto and Capella |^, who 
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Figure 19: Jjip suppression factor due to 
absorption on nucleons as obtained in the 
Glauber approximation with an absorption 
cross section of 7.3 ± 0.6 mb, in comparison 
with the NA38/50 data for the S + U and 
Pb + Pb systems (from 0). The suppres- 
sion in S + U is nicely explained, but the 
Glauber model fails in the Pb + Pb system. 
The excellent agreement for S + U appears 
to leave no room for comover absorption in 
this system. 



obtain acceptable fits with the parameter choices c^g^'' = 6.7(7.3) mb, <f = 0.55(1.0) mb, 
and critical area density of wounded nucleons of dN^"^/dy = 1.15(2.5) fm~^. 

What these two analyses demonstrate is that the present body of data on J/ ip production 
inp + A collisions leave very little room for additional suppression effects in S + U collisions, 
at least within the framework of the Glauber model. In order to obtain sufficiently strong 
suppression effects in Pb + Pb, a strongly nonlinear dependence of the exponent of ( plf ) 
on the comover density (here modeled as a threshold effect) is then required. Whereas 
Kharzeev et al.0 argue that this indicates a qualitatively new mechanism (such as color 



deconfinement), the authors of [32| do not reach this conclusion. 

The results for the hadronic absorption rate described above arguably point toward a fail- 
ure of the Glauber model to describe the Pb + Pb data, at least if only hadronic interactions 
are considered.^ Whether the reason for this failure is the emergence of a new suppression 
mechanism, or whether the A-dependence of the effective comover density is strongly un- 
derestimated by the Glauber approach, is impossible to tell from the existing data. The 
latter possibility has been recently addressed in where it was pointed out that the co- 
movers responsible for J/ip suppression have to be "semihard" in order to resolve the color 
dipole structure of the J/ip. A perturbative QCD picture then suggests that the area den- 
sity of these comovers grows as (^1,^2)^^'^ in collisions at the SPS, much faster than the 
{A^^^ + A}^^^) growth predicted by the wounded nucleon model. Of course, such comovers 
would be better represented as partonic excitations, not hadronic ones. Indeed, a hybrid 

^Not all authors share this opinion. In particular, Gavin[Q and Capella[Q have maintained that ab- 
sorption by hadronic or "pre-hadronic" comovers can describe all observations. 
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parton cascade model predicts just such a dependence for the partonic components of 
secondary particles (see Fig. pO|). 
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Figure 20: Partonic energy densities for S + U and Pb + Pb in the parton cascade model 
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On the other hand, the suppression of the ip' state is so strong already in S + U collisions 
that its description requires comover suppression. This is almost trivially obvious from the 
fact that the ip' /{J/il)) ratio is constant m p + A collisions, but drops significantly and 
progressively with decreasing impact parameter (growing dET/dy) in S + U [0. This has 
been analyzed by many authors P, | 



|37| , ^ who all reach similar conclusions (see Fig. |21 
for an example). Values for the absorption of the ip' by comovers are obtained in the range 
(T^^ ('?/'') = 8 mb, which is quite reasonable in view of the rather large size of the ip' and of its 



low dissociation threshold. In fact, the authors of |^ argue that the ip' suppression observed 
in Pb + Pb is not quite as strong as expected on the basis of the S + U data, and that a 
regeneration mechanism must be introduced to obtain a good fit of the Pb + Pb data. 



4.2 Microscopic Models 

Microscopic models of heavy ion reactions allow us to study the question whether the models 
based on Glauber theory adequately describe the space-time dynamics of these reactions. 
This includes questions such as whether the time dependence of the comover density and 
the comover composition is modeled correctly in the Glauber approach and whether the 
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Figure 21: ?/^' suppression factor due to 
absorption on nucleons (with a cross sec- 
tion of 7.3 lb 0.6 mb) as well as absorption 
on comoving hadrons (with a cross section 
of about 10 mb), in comparison with the 
NA38/50 data for the S + U and Pb + Pb 
systems (from 0). Comover absorption is 
required in both systems if one wants to 
obtain agreement between the data and the 
Glauber model. 
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transverse expansion, which is neglected in the Glauber model, plays a significant role. 

An extensive study of such questions was recently made by Bratkovskaya and Gassing in 
the framework of a relativistic hadronic cascade model The model treats the nuclear 
collision as a sequence of binary hadronic scatterings with the possibility of a formation time 
for all newly created particles. When a J/ip particle is produced in a primary nucleon-nucleon 
collision, it can interact with another hadron and be dissociated. The authors assume that 
the (cc) pair is originally produced in the color octet state and converts to a color singlet 
J/ip after Tf = 0.7 fm/c. 

It is useful to begin taking a look at the general space-time structure of a nuclear collision 
event. The Pb + Pb system differs from higher systems (p + Pb, S + U) in that a significant 
fraction of J/ip are formed as color singlets during the time period in which the primary 
nucleon-nucleon collisions occur. This is so because at the GERN-SPS energy, corresponding 
to 7cm ~ 9, the Pb nuclei are Lorentz contracted to a longitudinal width of only about 1.5 
fm, exceeding the assumed J/ip formation length by a factor two. Accordingly, the J/ip 
states are on average produced in a more violent environment than in the case for lighter 
collision systems. The density of comovers is higher in Pb + Pb than in S -|- U and falls off 
more slowly. 

The Gassing-Bratkovskaya model predicts a peak density of pions which slightly exceeds 
1 fm^ in the Pb + Pb system. A look at our Fig. |15| shows that this pion density is equivalent 
to a temperature T ^ 300 MeV under equilibrium conditions. The same holds true for the 
peak density of p-mesons (about 0.6 fm~^), indicating that the comover abundances may be 
close to chemical, if not thermal, equilibrium. In other words, the model predicts hadronic 
comover densities far in excess of those in the range of the presumed validity of the hadronic 
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phase. The peak densities predicted for the S + U system are lower and correspond to those 
of a thermal environment of "only" T ^ 250 MeV. 
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Figure 22: Density of comovers as function 
of time in the cascade model of Gassing and 
Bratkovskaya. 
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Figure 23: Time distribution of J/i/j ab- 
sorption events. Most absorptions occur 
within 3 fm/c after the onset of the reac- 
tion. 



Figure |23| shows that most of the absorption of J/ip by comovers in this model occurs 
within the first 3 fm/c after the initial reaction. At that time (t = 5 fm/c in Fig. |23| ), 
the comover density still corresponds to a thermal bath at T ^ 250 MeV in Pb + Pb 
and T ^ 210 MeV in S + U. This confirms the conclusion reached in Section III.C, that 
significant dissociation rates for J/-?/; in a thermal hadronic environment require temperatures 
far beyond T^, at least in the range T ^ 250 — 300 MeV. 

A second interesting observation is that the time interval during which the J/ip can be 
efficiently dissociated is so short that transverse expansion effects can be safely neglected. 
Hence, the Glauber approximation should be an excellent approximation, and in many re- 
spects preferable to microscopic cascade models, because it allows for a more consistent and 
economical treatment of quantum effects.[| 

^Cascade models can be useful tools for estimating quantities such as the effective optical thickness of the 
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In conclusion, hadronic cascade models provide useful tools for testing the consistency 
of the hadronic comover suppression scenario. At present, their results confirm that this 
scenario has difficulty explaining an observed "anomalous" suppression within the range of 
validity of the hadronic gas model. 



5 Addendum: Developments since June 1997 

Since mid-1997, a large number of articles have appeared addressing the implications of the 
new data from experiment NA50 on J/ip and ip' production in Pb + Pb collisions. This 
activity was triggered by the new results presented by the NA50 collaboration, first at the 
Quark Matter 1997 conference in Tsukuba |]14| and later at the Moriond meeting in March 
1998. The new data from the analysis of the 1996 Pb beam run at CERN, reproduced in 
Figures || and ^ revealed a pronounced sudden drop in the ratio of J / ip production to Drell- 
Yan pairs in the Pb + Pb system at an impact parameter around 8 fm (corresponding to 
a nuclear absorption length L ^ 8 fm). Furthermore, the new data confirmed the presence 
of an "anomalous" suppression effect in Pb + Pb, and their much higher statistics leaves 
hardly any doubt that the effect is real. 

After a careful reanalysis of the viability of hadronic comover models Vogt P6[ concluded 



that the data for p + A and S + U collisions leave insufficient room for hadronic comovers to 
make an explanation of the suppression of J / ip in Pb + Pb in terms of hadronic absorption 
viable. This agrees with the results of Kharzeev, et al. §]. Vogt finds no difficulty with an 
interpretation of the ip' data from NA50 in terms of hadronic absorption. Vogt's analysis 
includes the suppression effects on the feed-down to J / ip from the excited states Xc (normally 
about 30%) and ip' (12%). The S -|- U data can be fitted with a'^^'' = 7.3 mb and no comover 
suppression or with a^,^^ = 4.8 mb and (Xco = 0.67 mb, 1.6 mb, and 2.5 mb for J/ip, Xc, 
and respectively. The extrapolation to Pb + Pb then fails to explain the NA50 data, 
if the comover density is scaled as Et, as expected on the basis of multiparticle production 
models, such as the wounded nucleon model. A reasonable agreement with the Pb + Pb 
data would require a much faster than linear rise of the comover density with Et, such as 

E^T^\f\ 



Conclusions similar to those of Vogt are reached by Wong E^, who analyzes in detail 



comoving matter, J drvaahsdN {t) / dy ^ as input into simplified Glauber calculations. However, one needs to 
exercise caution since quantum coherence effects that may affect N{t) are neglected. 

®Such a fast rise could be explained by assuming that the comover density relevant for J absorption 
scales like other hard QCD processes, as suggested by the arguments presented in pO|. 
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some implications of the competing color singlet and octet production mechanisms for the 
J/ip and . Martins and Blaschke [^], based on their hadronic charm exchange model |p7|| , 
also rule out the hadronic absorption scenario. 

On the other hand, Dias de Deus and Seixas argue that a very simple, schematic 
model of comover suppression fits the NA50 data in a roundabout way; however, their model 
does not explain the apparently different behavior of the low- and high-E'^ regions in the Pb 
+ Pb system. Frankel and Frati argue that the observed suppression can be explained in 
terms of the energy loss of partons (gluons) as they travel through the colliding nuclei. Apart 
from being unable to describe a sudden drop in the J / ip yield, it is not entirely clear how this 
model can escape the stringent limits on energy loss set by the Drell-Yan phenomenology. 

A number of authors have quantitatively investigated the effectiveness of various absorp- 
tion mechanisms for J/ip in dense matter. The reaction TT + J/ip^ip' + Trm the presence of 
a thermal pion gas was studied in two articles. Chen and Savage analyzed the reaction 
within the framework of chiral perturbation theory, finding a very small average cross section 
of order 0.01 mb at thermal energies. Sorge, Shuryak, and Zahed [^| made the assumption 
that the tctt spectrum in the decay ip' —>■ J/ipTTir is dominated by a scalar resonance. In 
their model the ratio ip' /{J/ip) can chemically equilibrate to the value 0.05 observed in the 
highest Et bins of the Pb + Pb and S + U data, if the mass of the scalar resonance drops at 
high density, enhancing the reaction rate. Shuryak and Teaney |Q calculated the reaction 
IT + J/ip ^ rjc + p, but found that the thermal rate was only about 10~^ (fm/c)~^. 

Several authors PB|, ESI, EM have emphasized that the major part, if not all, of the 



additional suppression of the J/ip yield seen in the Pb + Pb data can be understood as 
virtually complete absorption of the fraction (about 30%) of J/ip which originates from 
decays of Xc- The deconfinement scenario favors this explanation, because the Xc disappears 
almost immediately at T^, while the J/ip survives to higher temperatures. This picture 



is supported by recent, improved lattice calculations |56[. QCD-based calculations of the 



absorption cross sections of J/tp and Xc on nucleons also show the higher vulnerability of 



the Xc state at low energies ^7j. Gerland et al. |^ argue that the stronger absorption of 
Xc is already needed to explain the p + A data and hence must play a significant role in any 
explanation of the heavy ion data for J/ ip suppression. 

The picture of J/ip suppression by QCD strings has been embraced with considerable 



enthusiasm lately. Following [^, Geiss et al. |Q simulated this effect within the HSD 
model discussed in section 4.2. They found reasonable agreement with the NA50 data by 
setting the string radius to 0.2-0.3 fm and assuming that any J/ip immediately dissociates 
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when entering a string. With a string area density of about 3/fm^, this imphes that at least 
half the transverse cross section of the reaction volume in Pb + Pb is filled with QCD strings. 
A related, but more schematic study by Nardi and Satz concludes that the sudden drop 
in the J/ip yield observed by NA50 coincides with the point where QCD strings become 
so dense that large overlapping clusters develop for an assumed string radius of 0.2 fm. In 
this picture the drop is associated with a color percolation phase transition which causes the 
disappearance of the Xc state. The existence of such a phase transition was also emphasized 



by Braun, Pajares, and Ranft ||60|, who studied the percolation of strings analytically as well 
as by simulations. Much earlier studies of J/ip absorption reached similar conclusions [pO|] . 



6 Summary and Outlook 

Have heavy quark bound states, the J/ip and ip' , fulfilled their promise as probes of the 
structure of dense hadronic matter formed in relativistic heavy ion collisions? Although the 
final answer to this question is still extant, several conclusions can be drawn now. Clearly, 
the original expectation that any significant reduction of J/ formation in nuclear collisions 
would signal the creation of a quark-gluon plasma had to be revised. We now understand 
that the suppression observed with light nuclear projectiles at the CERN-SPS is a smooth 
extrapolation of the less-than-linear target mass dependence found in p + A reactions. We 
also know a credible mechanism for this suppression effect: the color-octet formation model. 
However, it needs to be stressed that we have almost no direct evidence for the color-octet 
(cc) state, and we lack a quantitative description of the transition of the color-octet pair to a 
color-singlet charmonium state from first principles. For these reasons, the color-octet model 
must be considered a reasonable, even likely, explanation of the suppression from p + Ato S 
+ U reactions, but it is far from being firmly and quantitatively established. 

The available experimental evidence, combined with theoretical arguments, points toward 
the emergence of novel J/ip suppression mechanism in Pb -|- Pb collisions. Because no 
experiments with projectiles between ^^S and ^osp]^ have been done, it is unclear whether 
this new mechanism sets in suddenly or gradually. However, the magnitude of the effect is 
such that it can hardly be explained as the absorption of color-singlet J/ mesons on hadronic 
comovers, unless one wants to invoke the existence of hadrons, as we know them, at densities 
exceeding 1 fm~^. Moreover, a sizable comover suppression effect should be visible already 
in S + U collisions, if the "anomalous" suppression seen in the Pb + Pb system could be 
accounted for by absorption on hadrons. On the basis of the existing experimental evidence. 
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the structure of comoving matter has to undergo a significnt change between S + U and Pb 
+ Pb, if comover absorption is at the origin of the suppression effect in Pb + Pb coUisions. 

In the following, we enumerate some weaknesses of the present arguments regarding the 
"anomalous" J/ip suppression as evidence for the creation of a quark-gluon plasma in Pb 
+ Pb collisions at CERN. We also list some issues warranting further experimental and 
theoretical study. 

1. The argument that J/ip suppression in S + U collisions is "normal" but that seen in 
Pb + Pb collisions has a novel, "abnormal" component relies critically on the precision 
of the extrapolation of the p + A results. As the advocates of hadronic comover models 
have pointed out, the Pb + Pb data would lose their cogency if the new suppression 
mechanism were found to set in gradually as the nuclear projectiles get heavier. How 
large is the error in the value cr^g^^'* = 7.3 mb, i.e. how uncertain is the extrapolation of 
the Gerschel-Hiifner line to nuclear systems? The extrapolation from p + A to nuclear 
collisions depends critically on the assumption that the time for color neutralization is 
sufficiently long, so that the (cc) pair remains in the color octet state while the nuclear 
matter sweeps by. Can one really exclude cr^,^^ = 5 mb if all errors in the p + A data 
and the extrapolation to nuclear systems are taken into account? 

2. The viability of the hadronic comover model critically depends on the value of the 
absorption cross section for J/ ip on hadrons. It is possible to obtain a safe upper bound 
for cr^l^^'^ in a hadronic gas of imprecisely known composition? A systematic study of 
J/ip absorption within effective meson theories including mesons containing heavy 
quarks would be extremely useful. Another important study would include heavier 
hadronic states in the operator-product expansion apaproach to J/ ip absorption. Since 
{h\asE'^\h) ~ ml, heavier mesons can act as ample sources of gluons that may excite 
and dissociate the J/ijj. (The higher dissociation efficiency of p-mesons is also visible in 
the D-meson exchange model, see section III.C.) For example, it would be interesting 
to explore the J/ ip dissociation rate in a resonance gas with a Hagedorn-type excitation 
spectrum. 

3. An important issue that has not been investigated extensively is the problem of Xc 
absorption [^. Clearly, must lie somewhere between the cross sections for J/ip 
and ip' absorption. The question is if and where "anomalous" Xc absorption sets in. 
Can hadronic comover absorption of Xc in S + U be ruled out? Can the "anomalous" 
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absorption of J/ip in Pb + Pb be attributed entirely to simiarly "anomalous" absorp- 
tion of the component due to feed-down from Xc? Vogt's extensive study |^6[ seems to 
indicate that it is difficult to find a scenario including feed-down from Xc to J/ip that 
is compatible with all existing heavy ion data. 

4. There are many improvements that could be made to microscopic models of J/ip sup- 
pression. Including a realistic energy dependence of (T^^'lf', similar to the one found in 
section III.C, would probably lead to an effective nonlinear dependence of the comover 
effect on dN/dy, because higher comover densities are usually correlated with increased 
average kinetic energy per particle. 

Because the p-r-spectrum of J/ ^/'-suppression will be used as a tool for discriminating 
between different mechanisms, a study of final state changes in the pr-distribution of 
J/ip, e.g. due to elastic scattering on comovers, would be interesting. In more general 
terms, microscopic transport models could be used to systematically assess kinematical 
effects that can invalidate the Glauber approximation in the context of J/ ip production, 
in particular, effects of transverse expansion and changes in the pt distribution of the 



J/^p by elastic scattering in the medium Both effects increase for heavier nuclei 



Desirable improvements in the experiment include a reduction of statistical and sys- 
tematic errors by taking more data. Higher statistics data would also permit the 
determination of the J ftp suppression factor in a variety of Et bins and over a range 
of values of px, the transverse momentum of J/ip- A systematic study of the mass 
dependence of the effect for several symmetric systems from S + S to Pb + Pb would 
be especially important, so that the specific predictions of the deconfinement model 
can be checked against the data. An excitation function would also be of interest, 
but due to the steep beam-energy dependence of J/ ip production it may be difficult to 
eliminate systematic errors with the required accuracy. As pointed out by Kharzeev 



and Satz |57], the inverse kinematic experiment (Pb + S) would allow for a test of 



formation time effects on the suppression factor. 

6. The comparison of the experimental data with theoretical models will be facilitated, 
if the data are plotted against measured variables (such as Et) rather than model 
dependent derived quantities (such as L). Of course, a detailed understanding of the 
experimental trigger conditions is indispensible for a meaningful comparison between 
experiment and theory. 
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A clear set of ground rules needs to be established for the data analysis as well as 
for the comparison with theory. If the ratio between J/il^ yield and Drell-Yan (DY) 
background is plotted, only those mass ranges should be included where the DY process 
can be unambiguously identified. It makes little sense to include the region under the 
J/%1) itself, where the DY pairs make only a small contribution. 

Theorists also need to take a careful look at the soundness of their models. While 
microscopic transport models may have the appearance of great realism, their quan- 
titative results can depend on many implicit model assumptions that are difficult to 
analyze. Thus, more schematic models will continue to provide useful analytic tools. 
However, blind or ad-hoc parameter fitting can easily confuse the situation. Parameter 
values must be related to other experimental data, providing support for the fit, and 
should not exceed the limits of applicabihty of the schematic model itself. In some 
cases, it may be possible to check by ah initio QCD calculations whether a certain 
parameter set makes sense. As the picture of J/V' suppression in nuclear reactions 
takes on a firmer shape, models must increasingly be judged by their ability to predict 
new observables quantitatively. 

Finally, with the experimental emphasis shifting from the CERN-SPS to RHIC in the 
years ahead, it may be possible to observe not only the nuclear suppression of the J/V', but 
also the suppression of the T in heavy collision systems. Parton cascade models predict initial 
energy densities of the order of 50 GeV/fm^, far exceeding the deconfinement threshold of the 
{hb) ground state. The advantage of the T would be that its suppression in p -|- ^4 collisions 
is much weaker than that of the J/tl) and, therefore, its almost complete suppression in a 
deconfined phase would be quite spectacular. 
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Appendix A: J/^ Dissociation in a Medium 



Here we consider the dissociation of a J/ip state in a medium within the Bhanot-Peskin 
formahsm. The fundamental idea is that the dissociation is initiated by the absorption of a 
gluon from the medium which excites the color-singlet (cc) pair into a color-octet continuum 
state \{cc)^^\e), where e denotes the final state energy. The S-matrix element is 

% = ^ / dt{ic-cf^e\gf-E{r,t)\iccY'^). (26) 

til J~oo 

We will make use of the dipole approximation assuming that the dissociation is dominated 
by long-range gluonic interactions. In the heavy quark limit, the binding energy of the (cc)'-^-' 
state is |a^mQ, whereas its inverse size is = la^mg. When -C 1, the characteristic 
wavelength of a gluon that is energetically capable of dissociating the J/tp state is larger 
than the J/ ip radius. This approximation is somewhat marginal for the J/ ip state; it should 
be very good for the T. After performing the SU(3)-color algebra, the dipole matrix element 
is given by 

-^|rfVy,^(r)rXr), (27) 

where ^piir) and ^%{r) are the spatial wavefunctions of the singlet and octet states, respec- 
tively. We neglect recoil effects, and we use a Coulombic Is-wavefunction for the singlet 
state and a plane wave for the octet state. The resulting S-matrix then takes the form 




2 ■ P 



^'■"TiV-^^^ITW- '''' 

Here E"-{uj) is the spectral distribution of gluon fields in the medium at frequency uj = eg — ei, 
a is the Bohr radius of the (cc)*^^-* state, pis the relative momentum of the (cc)'-^^ pair, and 
V denotes the quantization volume. Integrating over p and assuming color neutrality of the 
medium, i.e. 



{Et{u;)E'^{u;)) = -6,,6,,mu;)\'), (29) 

where i,j denote spatial vector and a,b color indices, we finally obtain the transition prob- 
abihty 

P,f = ^7ia,a' {\E{u;)\'). (30) 

Since (|i?(cij)p) is proportional to the total interaction time Tjj, the result (^) corresponds 
to a constant dissociation rate. 

The color-electric power density 7^(\E(u)\'^) of the medium can be evaluated analytically 
for a variety of simple models. One such model is a dilute gas of color charges with a given 
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density and momentum distribution. Denoting the Casimir of the color charge by Q^, such 
a model yields 

-L(|E(,.)|2)^^«.gV», (31) 

where p{oj) is a w-dependent, weighted average of the density of charges in the medium. 

If the medium is represented as a dilute gas of hadrons, the power spectrum is related 
to the gluon structure functions of the hadrons Gh{x) where x = w/ph, in the spirit of the 
calculation of Kharzeev and Satz [Q. Given a momentum spectrum of hadrons f{ph) and 
chemical abundances of different hadrons, the power spectrum can be evaluated. 

Crude estimates of the dissociation rate in these two models yield rather long survival 
times of the J/ip in any realistic hadronic medium. For example, counting the valence quarks 
in mesons and combining (^) with (pT]), would yield a dissociation rate 

Tdis ~ -TT^OlWPmeson, (32) 

y 

where only sufficiently energetic mesons are counted. 



Appendix B: Matrix Elements for Charm Exchange 

In this Appendix we list the explicit expressions for the invariant matrix elements for the 
charm exchange reactions described by the Feynman diagrams in Fig. |T^. We begin with 
the two diagrams for J/ip absorption on pions (diagrams |T^,b). Averaging over initial and 
summing over final isospin and spin states, we obtain 

^2 _ ( ^2 _ iP^yPnl 

1 '-'2 2 V /V '"vi / 

\Ma\^ = -g^DDdnDD* (g2 _ rn\)Y ^^"^^ 

where q = p-^ — Pd* = Pd ~ Pip is the momentum transfer of the reaction. Introducing the 
Mandelstam variables s, t and u, and the shorthand notation t' = t — m|,, we obtain 

= \91dd91dd' i U< - ml, + 2t' - (Ami '^1 + 2^' - ^\ ■ (34) 

The result for the diagram (|l^b) is identical, i.e. 



|Mfe|2 = |M,p. (35) 

The contributions of the diagrams (0c, d) need to be added coherently, because they lead 
to identical final states. However, their results are related by crossing symmetry {t ^ u). 
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After somewhat lengthy algebra one finds: 



uV \ " "ll I \ " ml ) 

+ ((' « u')] (36) 

where again t' = t — m'jj, u' = u — m'jj. Of course, the crossing related diagrams ([T3|c,d) yield 
the same contribution to the total absorption cross section, with an additional interference 
term. 

Finally, for completeness, we give the expression for the decay width of D* Dir: 
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